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Abstract 

Background: Adaptive radiation is the process by which a single ancestral species diversifies into many 
descendants adapted to exploit a wide range of habitats. The appearance of ecological opportunities, or the 
colonisation or adaptation to novel ecological resources, has been documented to promote adaptive radiation in 
many classic examples. Mutualistic interactions allow species to access resources untapped by competitors, but 
evidence shows that the effect of mutualism on species diversification can greatly vary among mutualistic systems. 
Here, we test whether the development of obligate mutualism with sea anemones allowed the clownfishes to 
radiate adaptively across the Indian and western Pacific oceans reef habitats. 

Results: We show that clownfishes morphological characters are linked with ecological niches associated with the 
sea anemones. This pattern is consistent with the ecological speciation hypothesis. Furthermore, the clownfishes 
show an increase in the rate of species diversification as well as rate of morphological evolution compared to their 
closest relatives without anemone mutualistic associations. 

Conclusions: The effect of mutualism on species diversification has only been studied in a limited number of 
groups. We present a case of adaptive radiation where mutualistic interaction is the likely key innovation, providing 
new insights into the mechanisms involved in the buildup of biodiversity. Due to a lack of barriers to dispersal, 
ecological speciation is rare in marine environments. Particular life-history characteristics of clownfishes likely 
reinforced reproductive isolation between populations, allowing rapid species diversification. 

Keywords: Ecological speciation. Diversification, Comparative method. Evolutionary rate, Brownian Motion, 
Pomacentridae 



Background 

The concept of adaptive radiation has been central to 
evolutionary biology since Darwin's work on Galapagos 
finches [1-3]. The general understanding of this process 
is that rates of ecomorphological changes and species di- 
versification will be increased by ecological opportunities 
offering available resources untapped by competing spe- 
cies [4]. Ecological opportunity can arise for four main 
reasons [5], the most widely described being the colon- 
isation of geographically isolated areas with depauperate 
fauna (e.g. cichlid fishes in East- African Great Lakes [6]). 
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The process is similar in the aftermath of a mass ex- 
tinction event, which allows surviving species to radiate 
rapidly by filling the available empty niches [7] . Modifi- 
cation of a resource can also trigger native species radi- 
ation as demonstrated by the radiation of Lupinus in 
high-elevation habitats that appeared during the An- 
dean uplift [8]. Finally, the appearance of a trait allow- 
ing new interactions with the environment, or key 
innovation, can create an opportunity for species radi- 
ation [9]. For example, the evolution of antifreeze 
glycoproteins found in notothenioid fishes of Antarctica 
is thought to have triggered their adaptive radiation by 
allowing survival in extreme environments [10]. In an 
analogous manner to key innovations, the evolution of 
mutualistic interactions can provide access to previ- 
ously inaccessible resources. For instance, phytopha- 
gous insects host mutualistic microbes, which enable 
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the breakdown and digestion of plant compounds by 
the insects [11]. While a plethora of case studies show- 
ing adaptive radiation driven by ecological opportunity 
offered by one of the aforementioned possibilities exist 
[9], examples involving mutualism are scarce (but see 
[12]). Since hosts shifts have allowed ecological speci- 
ation in a wide range of organisms, including coral- 
dwelling fish [13,14], there is a possibility for ecological 
speciation to occur in mutualistic systems. However, 
results from empirical and theoretic studies give contra- 
dictory evidence on the effect of mutualism on species 
diversification [15,16]. The topic is thus still debated 
and in need of further case studies. 

The clownfishes (or anemonefishes; subfamily Amphi- 
prioninae) are a group of 30 species within the damselfish 
family (Teleostei; Perciformes; Pomacentridae; [17]) and 
are emblematic species of coral reefs (Figure lA & IC). 
Their distribution spans from the Indian to the western 
Pacific Oceans (Figure IB) with their highest species rich- 
ness found in the Indo-Malay archipelago where up to 
nine species have been observed in sympatry [18]. Their 
complex association with sea anemones is now a textbook 
example for mutualistic interactions [19-21]. Clownfishes 
are left unharmed by the otherwise lethal nematocysts of 
the sea anemone tentacles. This ability is thought to come 
from a protective mucous coat that prevents the discharge 
of the nematocysts [22] and allows clownfishes to settle in 
sea anemones. The protection against predators pro- 
vided by the sea anemones is a direct advantage for 
clownfishes. Likewise, clownfishes chase the predators 



of the sea anemones. Furthermore, waste ammonia 
excreted by the clownfishes is used by the endosym- 
biotic dinoflagellates living in the sea anemone tis- 
sues, which makes it a three-way interaction [19,21]. 
The efficiency of the protection provided by the sea 
anemone is demonstrated by the extraordinary life span 
of clownfishes (ca. 30 years recorded for Amphiprion 
percula), which is twice as long as any other damselfish 
and six times greater than the expected longevity for a 
fish of that size [23]. 

While species of clownfishes can develop mutualistic 
interactions with up to ten species of sea anemones 
(Table 1), a large variation in host usage exists within the 
clade [20]. Eight host sea anemones have a widespread dis- 
tribution and two {Heteractis malu and Macrodactyla 
doreensis) have more restricted ranges, but are distributed 
around the centre of diversity for the clownfishes, making 
interaction between most clownfishes and host species 
geographically possible [20]. Although geographically 
widespread, sea anemone species differ in their preferred 
habitat (e.g. reef zonation, substrate, depth; [24]). It was 
shown that coexistence of multiple clownfish species was 
possible because of difference in host and habitat utilisa- 
tion [18]. It is therefore possible that the appearance of 
mutualism was the key innovation that allowed the clown- 
fishes to diversify in ecological niches associated with the 
different sea anemones species. However, this hypothesis 
has never been tested thoroughly. 

Examples of ecological speciation events are rare in 
marine ecosystems. This is likely due to the fact that 




Figure 1 Clownfishes and sea anemones mutualism, and geographic distribution. Illustration of the mutualistic relationship between 
Amphiprion chrysopterus and Heteractis crispo (A). The distribution of the damselfishes in blue and of the clownfishes in orange is shown in panel 
B. As for every clownfish species, the female Amphiprion percula (on top of the picture of panel C, here with Stichodactyla gigonteo) is bigger 
than the male beneath. 

) 



Table 1 Interaction matrix between clownflshes and their sea anemone hosts 



Cryptodendrum Entacmaea Heteractis Heteractis. Heteractis Macrodactyla Stichodactyla Stichodactyla Stichodactyla Heteractis 
adhaesivum quadricolor aurora crispa magnifica doreensis gigantea haddoni mertensii malu 
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Known interactions are shown by plus signs. Field records are lacking for A chagonsensis and A. thiellei, the most probable host is shown between parentheses. The species status of A leucokranos and A. thiellei is 
debated as they may be natural hybrids [20]. 
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long-distance dispersal is common among marine or- 
ganism [25]. However, clownfishes are known to have 
high larval retention to natal reefs [26,27]. They are 
also known to produce species-specific calls [28,29] 
that differ among geographic populations [30]. Such 
properties are likely to have reinforced reproductive 
isolation by reducing gene-flow between clownfishes 
populations and thus facilitated ecological speciation 
processes in clownfishes. 

In this study, we test whether the evolution of the mu- 
tualism with sea anemone in the clownfishes lineage is a 
key innovation that led to ecological adaptive radiation 
(sensu [4]). We ensure that the clownfishes are monophy- 
letic by building a phylogeny for the Pomacentridae family. 
Next, we show the occurrence of rapid speciation in the 
clownfishes by testing whether their diversification rate is 
higher than that of the other damselfishes. We then use an 
ordination method on the mutualistic interactions to de- 
scribe potential ecological niches associated with the sea 
anemones. We further apply phylogenetic comparative 
methods to test the association between morphological 
traits and the putative ecological niches. We finally measure 
the rate of evolution of the morphological traits to see if 
they fit the theoretical expectation of faster morphological 
evolutionary rate [4,31]. 

Methods 

Mutualism and clownfishes phenotype 

We collected data on the distribution of clownfishes among 
the 10 possible sea anemone host species (Table 1; 
[19,20,32]). We applied a multiple correspondence analysis 
(MCA) on the matrix of mutualistic interactions between 
sea anemones and clownfishes. The MCA analysis is the 
counterpart of principal component analysis for categorical 
data, which shows the underlying structure present in the 
dataset. We used the axes of the MCA that explained most 
of the variance in the subsequent analysis. This allowed us 
to determine in a multivariate space, the characteristics of 
the ecological niches used by the clownfishes and provided 
by mutualistic interactions. 

We extracted morphological measurements of the 
damselfish from the literature (mainly from [33], other 
sources are listed in the Additional file 1). It is thought 
that adaptation to sea anemones required modifications 
of the general shape as well as a change in swimming 
ability in host specialised clownfish. Indeed, continuous 
and fast swimming is not needed anymore because specia- 
lised species never venture far from their host [34]. We 
thus collected traits in the literature for all Pomacentridae 
species present in our phylogeny that are linked with body 
shape and swimming abilities as well as trophic niche, 
which is generally linked with habitat in Pomacentridae 
[35]. This analysis resulted in a matrix of eight morpho- 
logical traits (maximum standard length, the ratio between 



standard length and the greatest body depth or "body 
ratio", the count of hard and soft dorsal-fin rays, the 
count of soft anal-fin rays, the count of pectoral-fin 
rays, the number of gill rakers present on the first gill 
arch and the number of scales which possess a sen- 
sory tube or "lateral-line scales"). Standard length and 
body ratio describe the overall fish shape, which has 
been shown to be linked with adaptation towards 
habitats with differing water velocity regimes [36,37]. 
Fin morphology directly influences fish locomotory 
ability [38] and gill rakers are used as a proxy for the 
differentiation along the pelagic-benthic trophic re- 
source axis [39]. The number of lateral-line scales is 
one of the more pronounced morphological differ- 
ences between the clownfishes and other damselfishes 
[40], and may be of importance in the ecological adaptive 
radiation. It was not possible to take into account intra- 
specific variation in our analysis and we recorded a single 
value per trait estimated as the mean of the values obtained 
from the literature. To diminish potential allometric effects, 
all traits were log transformed before further analysis. 

Phylogeny and divergence time estimation 

We assembled DNA sequence data for 196 Pomacentridae 
species (170/356 damselfishes, 26/30 clownfishes) spanning 
all genera in the family (Accession numbers available in 
Additional file 2). Three cichlid species {Aequidens 
rivulatus, Thorichthys meekU Tomocichla sieboldii) were 
included as outgroups [41,42]. The concatenated sequence 
matrix was 6945 bp long and composed of six mitochon- 
drial and three nuclear gene regions (12S, 16S, ATP6-8, 
COI, cytochrome b, ND3, BMP-4, RAGl & RAG2). 
Each DNA region was aligned separately with MUSCLE 
[43] and ambiguously aligned nucleotides were removed 
using Gblocks [44]. 

After visually checking the alignment, we used BEAST 
[45] to simultaneously infer the phylogeny and estimate 
divergence times. We used a relaxed clock model, draw- 
ing substitution rates from a lognormal distribution. We 
partitioned the alignment by gene as it outperformed an 
unpartitioned analysis in Bayes factors in a similar data- 
set [41]. We selected, using Akaike information criterion 
values (AIC), the substitution model that fits best 
each partition with the function "phymltest" available 
in the Ape package [46] in R [47] (see the model 
choice in Additional file 3). We used the only fossil 
calibration point available for the basal node of the 
Pomacentridae to obtain absolute divergence time 
estimates. The fossil that is the earliest record of 
Pomacentridae (Monte Bolca, Italy) dates back to 50 
million years (MY) [48], which we used as minimum 
age with a lognormal prior (mean = 2; sd = 1.2; prior 
5-95% = 51.03-103.2) following [41]. We selected a lognor- 
mal prior to allow the basal node of the Pomacentridae to 
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reach back to -105 MY, which is the probable age of the 
Perciformes [49,50]. We performed two parallel BEAST 
runs, each 5*10^ generations long and sampled posterior 
distributions every 1,000 generations. We checked the 
convergence of the two chains, optimal sampling of model 
parameters and estimated the burn-in length in Tracer 
[45]. After the removal of 10,000 trees as burn-in, we 
merged both runs and inferred a maximum credibility 
phylogeny using TreeAnnotator [45]. Finally, we resampled 
from the posterior distribution 100 trees to be used in fur- 
ther analysis. These time calibrated trees are hereafter re- 
ferred to as the distribution of chronograms. As BEAST 
also outputs phylograms having branch lengths given in 
expected number of substitution per site, we applied the 
same resampling procedure to get a random sample 
of 100 phylograms. This allowed us to test our hypoth- 
esis on two sets of phylogenies instead of possibly bias- 
ing our results by choosing arbitrarily a specific branch 
length unit [51]. 

Diversification rate 

We used the package Diversitree [52] in R to test whether 
mutualism with sea anemones is linked with an increased 
diversification rate in the clownfishes as would be expected 
under the key innovation hypothesis. We applied the BiSSE 
method [53], which evaluates jointly the evolution of a bin- 
ary character (here presence or absence of mutualism with 
sea anemones), speciation and extinction rates. As we do 
not have a complete sampling of the Pomacentridae, we 
used an extension of the method that deals with incom- 
pletely sampled phylogenies [52]. A one-rate birth-death 
model is fitted to the whole tree and compared, using AIC 
and Likelihood ratio test, with an alternative model allow- 
ing two separate rates of speciation and extinction for 
clownfishes and damselfishes species. In this particular case, 
the clownfishes are a monophyletic group nested within the 
Pomacentridae phylogeny. No known clownfishes species 
has lost the mutualistic behaviour and we therefore forced 
the loss of mutualism in the model (parameter qlO) to a 
fixed null value. We optimised the other parameters of the 
model (rates of speciation, extinction and probability of 
character change) by Maximum Likelihood estimation in- 
dependently on each of the 100 randomly sampled chrono- 
grams to account for phylogenetic uncertainty. The rate of 
diversification was calculated by subtracting the extinction 
rate from the speciation rate. 

Phylogenetic signal and phenotype-environment 
correlation 

We estimated the phylogenetic signal in the morphological 
data on each of the 100 phylograms and chronograms with 
the K [54] and \ [55] indexes as implemented in the 
Phytools package [56] in R. Assessing the phylogenetic sig- 
nal of a trait on both phylograms and chronograms can 



help choose which branch length unit will be the most ap- 
propriate for comparative analysis [51]. For both indexes, 
a value close to 0 is diagnostic of a weak or nonexistent 
phylogenetic structure, while values close to one are 
expected if the data follows a Brownian motion (BM) 
model of character evolution. We performed randomisa- 
tion tests for the K and a likelihood ratio test for X to 
test for an observed phylogenetic signal significantly 
greater than 0. We repeated the analysis by taking 
into account only the clownfishes and this time also 
measuring the phylogenetic signal of the four first 
axes of the host usage MCA. 

Past competition creating character displacement be- 
tween related species will result in phenotypes that are 
correlated with resource usage [4] . We assessed whether 
the morphological traits collected are linked to host 
usage in the clownfish by measuring the correlation 
between each of the first four axes of the MCA and the 
eight morphological traits. We used phylogenetic gener- 
alised least squares (pGLS) as implemented in the caper 
package in R [57]. The \ parameter, which models the 
phylogenetic dependency of species trait values [55] 
was estimated by Maximum Likelihood and the model 
was replicated over each tree present in the samples 
of phylograms and chronograms. We assessed if the 
morphological variables explained a significant part of 
the variance in the model by running an ANOVA on 
the pGLS output. 

Morphological evolutionary rate 

We measured the differences in rate of morphological 
evolution between clownfishes and damselfishes by com- 
paring the fit of a single rate BM model to that of a mul- 
tiple rate model. It has been shown that other models 
could better fit the data than BM especially in adaptive 
radiations [58,59]. We choose to use BM because our 
goal is solely to compare the relative rate of evolution 
between groups and not the actual trait values. The sin- 
gle rate model assumes that all lineages accumulate the 
same amount of morphological variance per unit of 
time while the multiple model allows clownfishes to 
have a different rate of evolution than the damselfishes. 
Both models were specified in the Phytools package [56] 
that implements the non-censored version of a typical 
BROWNIE analysis [60]. The best fitting model was 
selected according to sample size corrected AIC (AICc). 
We analysed each of the recorded morphological traits 
on the two sets of 100 trees randomly sampled from the 
posterior distributions of phylograms and chronograms. 

Results 

Phylogenetic inference and divergence time 

Our maximum credibility phylogenetic tree shows strong 
support for the monophyly of the clownfishes with a high 
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posterior probability (PP = 0.98) for the basal node of the 
clade (Figure 2). The monophyly of the clownfishes as well 
as the general tree topology that we recovered was congru- 
ent with previous phylogenetic trees of the Pomacentridae 
[17,41,61]. Most of the nodes of the tree were highly sup- 
ported (PP>0.95, Figure 2 and Additional file 4). 

The estimated age of the root node of the Pomacentridae 
(-55 MY) was similar to previous findings [41]. We in- 
ferred a root age of -19 MY for the clownfish clade and 
found the crown Amphiprion group, which holds most of 
the clownfishes species (25 out of 30 species), to have an 
age of ~7 MY (Figure 2).We did not include in this diver- 
gence time analysis the closure of the Isthmus of Panama. 
This calibration point, which constrains the split of the 
pairs Abudefduf concolor/taurus and Abudefduf troschelli/ 
saxatilis [41] was not used, because biogeographic informa- 
tion can be uncertain when constraining the age of species 
splits [62,63]. However, we recovered similar dates and con- 
fidence intervals as estimated in [41]. 

Speciation rate 

We tested whether the evolution of the obligate mutual- 
ism with sea anemones fitted the expectation of a key 
innovation and was linked with an increased rate of 
speciation in the clownfish. We observed that mutual- 
ism with sea anemones was linked with higher speci- 
ation, extinction and diversification rates (Figure 3). 
The model allowing distinct rates of speciation and 
extinction for clownfishes and damselfishes also explained 
the data significantly better than the simpler model where 
both groups have the same rate (median likelihood ratio 
test P = 0.02). The dating and phylogenetic uncertainty 
are taken into account in the final result (Figure 3) 
because we ran these analyses on a random sample of 
100 chronograms. 

Ecological niche of host usage 

The first four axes of the MCA explained 76% of the 
total variance in mutualistic interactions among clown- 
fish species and were kept for the following analysis 
(Figures 4 and 5, see Additional files 5 and 6 for factorial 
maps). Using knowledge on sea anemone habitat and 
ecology [24], we could interpret the principal axes of the 
MCA. The first axis (35% of variance) showed a gra- 
dient of differing host usage by segregating generalists 
clownfishes (positive values) that have interactions with 
many sea anemone species from specialists (negative 
values), which have a small range of possible sea 
anemone hosts. The remaining axes showed gradients 
Unked with habitat utilisation. Indeed, the second axis 
(15% of variance) separates clownfishes species inter- 
acting with sea anemones that live on different types 
of substrate (e.g. Heteractis aurora on sand and Entacmaea 
quadricolor on rock). The third axis (14% of variance) 



shows principally a depth gradient and the fourth axis 
(12% of variance) exhibits a gradient between sand 
dwelling sea anemone species living either among or 
away from coral reefs. Although this has not been formally 
tested, species that are close in the MCA (Figures 4 and 5) 
and thus similar in host usage, seldom co-occur in the wild. 

Phylogenetic signal 

In the whole Pomacentridae family, the morphological 
traits were generally highly conserved (K and \ close 
to 1, Table 2). On average, the phylogenetic signal was 
closer to 1 when measured on phylograms than on 
chronograms. This would suggest that, for comparative 
phylogenetic methods that apply the BM model of 
character evolution, using phylograms would give more 
accurate results [51]. We nevertheless ran all the subse- 
quent analysis on both kinds of phylogenetic trees be- 
cause there is only a slight difference in the measured 
phylogenetic signal of our data between the phylograms 
and chronograms. 

In contrast to the whole Pomacentridae phylogeny, the 
phylogenetic signal of the morphological traits and the 
first four axes of the host usage MCA measured only on 
the clownfish clade were relatively weak (K and \ close 
to 0; Table 2). Only one morphological character had K 
value significantly different from 0 (three for the X) while 
no MCA axis showed this pattern (Table 3). 

Phenotype-environment correlation 

We assessed the correlation between clownfish morpho- 
logical traits and putative ecological niches (as described by 
the MCA axes) with a pGLS. The results for the models 
were congruent between analysis of phylograms and 
chronograms (Table 4). The analyses using phylograms 
found that all traits but body ratio and lateral line 
scales had a significant relation with the first axis of the 
MCA. Only standard length was important when cor- 
related with the third axis of the MCA. Results were 
similar for chronograms, except for the pectoral-fin and 
hard dorsal-fin rays counts that did not significantly ex- 
plain variation in the first axis of the MCA (Table 4). 

Morphological evolutionary rate 

We measured, for each morphological trait, the rate 
parameter of a BM model of character evolution to 
assess if the appearance of mutualism was linked to 
an increased rate of morphological evolution in the 
clownfishes. We estimated, using AICc, if a model 
where clownfishes and damselfishes have distinct rates 
explains the observed data better than did a model of 
common rate between the two groups. We found that, 
on phylograms, all traits studied had a larger rate of 
evolution in the clownfishes than in the damselfishes 
(Figure 6). The pattern was more variable when measured 
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Figure 2 Pomacentridae maximum credibility ctironogram. Outgroup taxa are shown in black, damselfishes in blue and clownfishes in orange. 
Error bars on node show the dating confidence intervals, scale is in MY. Numbers above nodes indicate Bayesian posterior probabilities. 



Litsios et at. BMC Evolutionary Biology 201 2, 12:21 2 
http://www.bionnedcentral.conn/1471 -21 48/1 2/212 



Page 8 of 1 5 



CD 

DC ^ 



□ One Rate Model 

□ Damselfishes 

□ Clownfishes 



Extinction 



Diversification 



Speciation 

Figure 3 Diversification analysis. Rates of speciation, extinction and diversification measured on tine sample of 100 clironograms. Rates of 
damselfislies (mean diversification rate = 0.09) are sliown in blue and clownfishes (mean diversification rate = 0.14) in orange. 



on chronograms. In this case, only the standard length, 
body ratio, soft anal fin-rays and gill rakers had a faster 
evolutionary rate in the clownfishes. To verify that the ele- 
vated rates found in the clownfishes were not due to the 
relatively short branches of the clownfish clade we simu- 
lated a continuous trait under a BM model with a single 
rate on our phylogenies. The rates we recovered were not 
different between clownfishes and damselfishes. 



Discussion 

We found that clownfishes exhibit patterns that are 
likely diagnostic of an ecological adaptive radiation 
via ecological speciation [4]. Following the acquisition 
of specific ability to interact and live with sea ane- 
mones, clownfishes diversified into multiple ecological 
niches linked with both host (Figure 7) and habitat 
use. Morphological evolution accelerated and distinct 



clownfish species developed convergent phenotypes 
correlated to the host-associated ecological niches. 



Pomacentridae phylogeny and diversification 

Our inferred phylogeny was congruent with previous 
work [17,41] and showed with high posterior prob- 
abilities that clownfish are monophyletic within the 
Pomacentridae family (Figure 2). We used the BiSSE 
method [53] to test if the appearance of mutualism 
acted as a key innovation and if it is associated with 
an increase in speciation rate. The BiSSE method is a 
powerful way to detect shifts in diversification rate 
linked with a binary trait, but the change of state of 
the binary trait does not usually correspond to a sin- 
gle monophyletic group as in our case. Nevertheless, 
we chose the BiSSE method as it takes into account 
the uncertainty in dating the appearance of mutualism 



C\J 
< 

O 



eg 

I 



■ S. mertensii 

^ A. chrysogaster " H. magnified 
^ A. ocellaris □ ^ gtyrora 
A. akallopisos^^i^^ „ ^ ■ 

^ o S. gigantea 

A.leucokranos Q^. ,,,,,,,3 ° Crispa 

□ S. tiaddoni 

A. latifasciatus^^ al<indynos ° ^' ^^''^^'^^'^ 

A niarioes ^ ^""^ ° ^- quadricolor 

A. nigripes '^p^ ^<^A. percula >^ , , 

A. sandaracinos r • ■ /-/. malu 

(5 (3 ^A.chrysopterus , C. adhaesivum 

A.allardi^^ A. perideraion 


A. latezonatus vV J 

^"*^A. seba 

A. rubrocinctus^^ 

P. biaculeatus 

A. frenatus f j 

A. chagosensis'^fl^^ 

A. mcculloclii 

A. barberi 

A. ephippium Vi^ 
A. 


melanopus 
^^^A. polymnus 

A. c/ar/c/;^^ 


-2 C 


) 2 4 6 8 

MCA 1 
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Heteractis: H, Cryptodendrum: C. 
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on the phylogeny. Indeed, mutualism with sea ane- 
mones likely appeared in the common ancestor of the 
clownfish, but the stem branch of the group is ~13MY 
long (Figure 2). It is therefore important, when estimat- 
ing speciation and extinction rates, to take into account 
the uncertainty in the time estimates that correspond 
to the appearance of this behaviour. We also con- 
strained the rate of mutualism loss to be null in the 
BiSSE model This takes into account the fact that all 
clownfish are nowadays mutualistic, but it may slightly 
bias our analysis as it is not impossible that a clownfish 
species will eventually revert to a non-mutualistic state. 
However, we do not think that this parameter will in- 
fluence significantly our results. Other methods exist to 

Table 2 Phylogenetic signal in the damselfishes 



Phylograms Chronograms 





K 


A 


K 


A 


Standard 
lengtli 


0.407±0.11 


0.934±0.09 


0.358±0.10 


0.883±0.04 


Body ratio 


0.205±0.08 


0.954±0.07 


0.233±0.09 


0.921 ±0.02 


Dorsal fin 
soft rays 


0.467±0.14 


0.925±0.06 


0.461 ±0.1 4 


0.886±0.02 


Dorsal fin 
hard rays 


0.513±0.10 


0.914±0.07 


0.578±0.12 


0.859±0.02 


Anal fin 
soft rays 


0.1 61 ±0.07 


0.911 ±0.07 


0.186±0.08 


0.872±0.03 


Lateral line 
scales 


0.856±0.30 


1.026±0.06 


0.973±0.33 


0.962±0.02 


Pectoral fin 
rays 


0.349±0.12 


0.842±0.06 


0.329±0.12 


0.842±0.04 


Gill rakers 


0.237±0.09 


1.011 ±0.07 


0.294±0.12 


0.949±0.02 



Blomberg's K and Pagel's X statistics of phylogenetic signal and standard 
deviation for the damselflsh morphological dataset In the posterior 
distributions of phylograms and chronograms. For each trait, the statistic 
closest to one Is Indicated In bold. 



infer speciation rates on phylogenies that do not need 
prior hypotheses on the location of the shift in rate 
[64]. In a recent paper [41], the likely nodes of diversifi- 
cation rate shifts were inferred for four coral reef fish fam- 
ilies using relative cladogenesis test [65] and MEDUSA 
[64]. The study included the Pomacentridae family and 
they consistently found across methods that the clown- 
fishes experienced a significant rate increase. The rate shift 

Table 3 Phylogenetic signal of morphology and MCA 
axes 



Phylograms Chronograms 





K 


A 


K 


A 


Standard 
length 


0.186±0.06 


0.578±0.07 


0.186±0.08 


0.647±0.06 


Body ratio 


0.149±0.09 


0.868±0.16 


0.147±0.10 


0.829±0.08 * 


Dorsal fin 
liard rays 


0.201 ±0.04 


0.1 31 ±0.24 


0.179±0.05 


0.085±0.20 


Dorsal fin 
soft rays 


0.169±0.07 


0.175±0.30 


0.150±0.07 


0.1 41 ±0.26 


Anal fin 
soft rays 


0.11 ±0.09 


0.941 ±0.1 4** 


0.109±0.10 


0.863±0.03 


Lateral line 
scales 


0.245±0.10 


0.686±0.14 


0.216±0.09 


0.666±0.21 


Pectoral 
fin rays 


0.245±0.13 * 


0.977±0.15 


0.239±0.14 


0.89±0.08 


Gill rakers 


0.068±0.05 


0.721 ±0.11 * 


0.068±0.06 


0.738±0.05 


MCA 1 


0.097±0.06 


0.322±0.41 


0.088±0.05 


0.326±0.40 


MCA 2 


0.033±0.03 


0±0 


0.032±0.03 


0±0 


MCA 3 


0.108±0.05 


0.172±0.12 


0.106±0.05 


0.285±0.12 


MCA 4 


0.04±0.03 


0.021 ±0.1 5 


0.038±0.03 


0±0 



Blomberg's K and Pagel's X statistics of phylogenetic signal for the clownfish 
morphological dataset and host usage MCA In the posterior distributions of 
phylograms and chronograms. Phylogenetic signal significantly different than 
0 Is signified by asterisk = P-values <0.05, = P-values <0.01). Significance 
was assessed with a randomisation test for the K and likelihood ratio test for X. 
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Table 4 Correlation between morphological traits and 
MCA axes 

Phylograms Chronograms 

Morphological Coefficients Error Coefficients Error 
trait 



Table 4 Correlation between morphological traits and 
MCA axes (Continued) 



Standard 


1.498 


1.247 


1.56 


1.248 


length 


±0.035 


±0.015 


±-0.060 


±0.010 


Body ratio 


9.968 


3.393 


9.783 


3.413 




±0.191 


±0.032 


±0.670 


±0.018 


Dorsal fin 


7.273 


6.502 


7.820 


6.337 


hard rays 


±0.567 


±0.097 


±0.112 


±0.180 


Dorsal fin 


0.651 


6.551 


1.395 


6.395 


soft rays 


±0.609 


±0.073 


±0.400 


±0.321 


Anal fin 


9.599 


7.148 


9.327 


7.142 


soft ray 


±0.288 


±0.065 


±0.908 


±0.030 


Pectoral 


-2.244 


4.883 


-2.427 


4.864 


rays 


±0.161 


±0.044 


±1.107 


±0.054 


Lateral-line 


-10.471 


4.058 


-10.593 


3.949 


scales 


±0.154 


±0.083 


±0.200 


±0.018 


Gill rakers 


-0.620 


5.352 


-0.742 


5.316 




±0.151 


±0.067 


±1.517 


±0.199 



The table shows results and standard deviation of pGLS. Results in bold 
indicate variables explaining significant variation in the dependent variable as 
shown by the ANOVA on the pGLS output. Median adjusted of the models 
on phylograms, MCA 1 = 0.75, MCA 2 = 0.15, MCA 3 = 0.33, MCA 4 = 0.23, and 
on chronograms, MCA 1 = 0.65, MCA 2 = 0.16, MCA 3 = 0.36, MCA 4 = 0.22. 

was either placed at the origin of the clownfish or at 
the base of the Amphiprion crown group. The fact that 
a method which does not need a priori information 
on the location of the diversification rate shift recov- 
ered similar result as in our study confirms the 
strength of the diversification rate shift that occurred 
in the clownfish clade. 

It should be noted that the extinction rate also 
increases in the clownfishes (Figure 3). A possible ex- 
planation is that during the diversification process, some 
clownfish lineages did not leave any descendants as they 
were ecologically replaced by other more competitive 
clownfish species. Such events could have occurred during 
the long branch that is basal to the Amphiprion crown 
group. Also, it has been suggested that elevated extinction 
rates in reef associated fishes could be linked with a poten- 
tial refuge effect of the coral reefs in the aftermath of pro- 
longed extinction events [41]. Yet, to be able to recover 
extinction rate with high confidence and test such hypoth- 
eses, one would need clownfishes fossils [66], which are 
not available. However, when compared to the damsel- 
fishes, the diversification rate of the clownfishes was still 
higher (Figure 3), showing that the extinction rate was not 
sufficient to slow down diversification [41]. 

Effect of mutualism and host-associated niches on 
clownfishes evolution 

We measured the phylogenetic signal of each morpho- 
logical trait on the samples of chronograms and phylo- 
grams for the Pomacentridae (Tables 2 and 3). All 
traits showed a signal close to one (the expected 



Standard 


0.264 


1.785 


-0.093 


2.188 


length 


±0.321 


±0.115 


±0.755 


±0.408 


Body ratio 


0.039 


4.480 


-0.036 


5.618 




±1.765 


±0.468 


±1.435 


±1.423 


Dorsal fin 


8.617 


11.096 


3.661 


1 1 .074 


hard rays 


±1.941 


±0.715 


±4.792 


±2.168 


Dorsal fin 


-6.904 


12.666 


-9.241 


12.270 


soft rays 


±2.517 


±0.820 


±3.278 


±1.496 


Anal fin 


0.155 


9.593 


-2.716 


11.252 


soft ray 


±2.183 


±0.81 1 


±4.671 


±3.342 


Pectoral 


-7.830 


6.951 


-4.235 


8.226 


rays 


±3.415 


±1.247 


±5.400 


±2.055 


Lateral-line 


-4.353 


6.168 


-1.190 


7.030 


scales 


±1.593 


±0.349 


±4.424 


±-1.230 


Gill rakers 


16.496 


5.591 


15.797 


7.761 




±1.708 


±0.443 


±2.008 


±2.955 


Standard 


-0.920 


1.452 


-1.012 


1.473 


length 


±0.085 


±0.013 


±0.057 


±0.031 


Body ratio 


3.505 


3.950 


3.781 


4.029 




±0.143 


±0.054 


±0.459 


±0.024 


Dorsal fin 


3.751 


7.569 


2.811 


7.484 


hard rays 


±0.518 


±0.115 


±0.488 


±0.325 


Dorsal fin 


10.376 


7.626 


9.263 


7.554 


soft rays 


±0.431 


±0.095 


±0.831 


±0.511 


Anal fin 


-16.490 


8.321 


-16.398 


8.433 


soft ray 


±0.543 


±0.074 


±0.890 


±0.060 


Pectoral 


7.134±0.241 


5.685 


7.356 


5.742 


rays 




±0.060 


±0.747 


±0.003 


Lateral-line 


0.024±0.295 


4.724 


0.253 


4.663 


scales 




±0.080 


±1.082 


±0.078 


Gill rakers 


4.083±0.344 


6.23 


4.429 


6.274 






±-0.080 


±1.397 


±0.189 


Standard 


-2.421 


1.233 


-2.439 


1.227 


length 


±0.097 


±0.017 


±0.179 


±0.013 


Body ratio 


1.134 


3.356 


1.024 


3.355 




±0.157 


±0.072 


±0.206 


±0.010 


Dorsal fin 


-1.219 


6.430 


-1.458 


6.231 


hard rays 


±0.472 


±0.134 


±0.415 


±0.194 


Dorsal fin 


-1.975 


6.479 


-2.173 


6.288 


soft rays 


±0.435 


±0.131 


±0.112 


±0.336 


Anal fin 


-1.51 


7.069 


-1.429 


7.023 


soft ray 


±-0.480 


±0.126 


±1.638 


±0.014 


Pectoral 


-5.235 


4.829 


-5.158 


4.782 


rays 


±0.191 


±0.090 


±0.877 


±0.043 


Lateral-line 


-0.05 


4.013 


0.088 


3.883 


scales 


±-0.300 


±0.091 


±1.333 


±0.026 


Gill rakers 


-1.585 


5.293 


-1.566 


5.227 




±0.452 


±0.093 


±0.030 


±0.189 
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Figure 6 Rate of morphological evolution. Evolution of the morphological traits measured on the sample of 100 phylograms (A) and 
chronograms (B). Rates of clownfishes are in orange, damselfishes in blue and the one-rate model is shown in white. The P-values of model 
comparison by likelihood ratio test is signified by asterisk f = P-values <0.05, = P-values <0.01). The schematic position of each morphological 
trait is shown on the clownfish drawing. 



outcome of BM evolution). However, when assessed 
only on the clownfish clade, phylogenetic signal 
dropped and only pectoral rays count had a K signifi- 
cantly larger than 0. Furthermore, the phylogenetic 
signal of the host usage MCA axes were never signifi- 
cantly different from 0. While interpreting an evolu- 
tionary process directly from this result can be 
problematic [67], a low phylogenetic signal can be 
found in lineages that show convergent adaptive evo- 
lution [68], which is likely the case in the clownfishes. 

We hypothesised that, following the appearance of 
mutualism, clownfishes radiated in the niches associated 
with the sea anemones. We described the most import- 
ant axes of variation in mutualistic interactions with an 
MCA. The first axis depicted the generalist-specialist host 
usage gradient, but all three other axes showed gradients 
linked with the habitat preferences of the sea anemones. 
Indeed, clownfishes that interact with sea anemones 
species living in similar reef micro-habitats (i.e. substrate 



type, depth) cluster together in the analysis. This suggests 
that clownfish species are first distributed along a generalist 
to specialist axis, and then, specialist clownfishes interact 
only with sea anemone species living in a particular habitat 
type. This has been shown in a previous empirical study 
[18], where clownfish species coexisting in a reef were dis- 
tributed in different habitats. Ecological sorting of clownfish 
species along the different ecological gradients linked with 
their hosts is what is expected if resource competition, 
which is the main driver of adaptive radiation, acted on the 
evolutionary process [4,69]. 

We tested if the observed resource partitioning in 
different ecological niches, likely due to past competi- 
tion between ecologically similar species, resulted in mor- 
phological adaptation to resource use (i.e. host and habitat 
use in clownfishes). We sequentially fitted each MCA axis 
to a set of morphological traits taking into account the 
phylogenetic relationships between species. We found that 
an important part of the variation in the MCA axes could 
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Figure 7 Chronogram of the clownfishes radiation. Branch lengths are given in MY. The interacting sea anemone species are shown for each 
clownfish species. Sea anemone names abbreviations as in Figure 5. 



be explained by the morphological traits of our dataset for 
the first and third axes (Table 4). Increasing number of 
hosts (represented by increasing values of the MCA 1 axis) 
was correlated with a bigger size, more hard dorsal 
and soft anal rays, and more gill rakers, while size 
was positively correlated with deeper habitats (repre- 
sented by decreasing values of the MCA 3 axis). In clown- 
fishes, the phenotype-environment correlation relates to 
both host usage (generalists/specialist gradient) and habi- 
tat (substrate, depth). This contrasts with many exam- 
ples of adaptive radiation, where the resource axis 
has usually a single dimension representing habitat. There- 
fore, mutualism can be seen as a a key innovation that 
offered untapped habitat for colonisation, but also allowed 
diversification to happen on the host usage resource axes. 

The morphological traits studied are primarily used 
for taxonomic purposes but they can still give important 
functional information for the evolution of the clown- 
fishes. Size and fin traits are related to the locomotion 
ability in various types of water velocities, while gill 
rakers can be used as a proxy for the trophic level. The 
picture given by our analyses is that generalist 



clownfishes (e.g. A. clarkii) will likely eat more plank- 
tonic food (and thus have more gill rakers) and be better 
swimmers than specialists, which never leave the close 
vicinity of their sea anemone host. Clownfishes that 
interact with sea anemones occurring at deeper depth 
have also a bigger size, allowing for better locomotion in 
areas where the water velocity is likely to be higher com- 
pared to more shallow and sheltered zones [70]. While 
our results show correlations between phenotype and envir- 
onment, we do not test for trait utility. This would require 
a strict experimental setting that was out of scope for this 
paper. More studies are definitely needed to better describe 
the adaptive advantage that those traits may provide in the 
ecological context of the mutualistic interaction. 

Following an ecological opportunity, the rate of mor- 
phological evolution is hypothesised to be elevated in 
the traits that are functionally related to the ecological 
niches filled during the radiation process [4,71,72]. We 
tested this hypothesis on the eight morphological 
traits studied and found that, on phylograms, they all 
evolved at higher rates in clownfishes than in dam- 
selfishes (Figure 6). The picture is similar when rates 
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are measured on chronograms, although only four 
characters are evolving significantly faster in the clown- 
fishes than in the damselfishes in this case. Following the 
comparison of phylogenetic signal that was made between 
phylograms and chronograms, phylograms are assumed to 
give, in this case, more accurate results [51]. The fact that 
all traits evolve at a higher rate is congruent with the pGLS 
results, which shows all but two traits (lateral-line scales 
and body ratio) being significantly linked with host and 
habitat usage. It is probable that lateral-line scales and body 
ratio evolutionary rates were accelerated in clownfishes be- 
cause they are indirectly correlated to the ecological niche 
through another trait. Thus they would not be themselves 
correlated to host-usage/habitat but still show accelerated 
rates of evolution. A clear followup to this broad descrip- 
tion of clownfishes morphologies would be to extend our 
analysis and use a morphometric approach (e.g. [36]) to be 
able to give an accurate description of the different clown- 
fishes ecotypes. Further studies will then be needed to test, 
in an experimental framework, trait utility, which is one of 
the diagnostic criteria of adaptive radiation [4]. The chem- 
ical biology of the interaction between clownfishes and 
sea anemones is also far from being solved [22]. It is 
thus possible that unknown characteristics associated 
for example, with the clownfishes mucus could be 
linked with the variation in interaction between clown- 
fishes and specific sea anemones. 

We did not take into account the distribution of the 
species in this study. Geographical isolation, coupled 
with ecological differentiation, could also be at the ori- 
gin of the evolutionary pattern found here. While sev- 
eral clownfish species are local endemics that likely 
originated through vicariance events rather than ecological 
speciation, the majority of the species (17 out of 30) have 
overlapping distributions centred on the Indo-Malay 
archipelago. If reproductive isolation was solely due 
to geography, the latter species would have likely dis- 
appeared through hybridisation, which easily happens in 
captivity [73]. Moreover, clownfishes that have similar 
MCA values usually do not overlap in geographical 
distribution (e.g. A. latifasciatus, A. nigripes and A, 
sandaracinos in Figure 4), and sister species always 
differ in host usage as can be seen on Figure 7. Such 
a pattern could indicate that, in a given biogeographic 
region, only one species per ecological niche can sub- 
sist, but also that ecologically similar species evolved 
independently in geographically separated areas. This 
outUnes the need for a thorough biogeographic analysis 
that would help clarify the effect of geography on the 
evolution of the clownfish. 

Conclusion 

Our study shows that clownfishes likely experienced an 
adaptive radiation through ecological speciation. The 



obligate mutuaUsm with sea anemones is thought to be 
the key innovation that allowed clownfishes to radiate 
rapidly in untapped ecological niches. As expected under 
the ecological theory of adaptive radiation [4], it increased 
diversification as well as rates of morphological evolution. 
Clownfishes experienced rapid and convergent morpho- 
logical changes that were correlated with the different eco- 
logical niches offered by the host anemones. In marine 
environments, barriers to dispersal are uncommon, 
which makes ecological speciation less likely than in 
more isolated landscapes [25]. However clownfishes 
show a very short dispersal period compared to other 
damselfishes [74]. In conjunction with a high reten- 
tion of larva to natal reef [27] and population specific 
calls [30], restricted dispersal likely reinforced repro- 
ductive isolation between clownfish species allowing 
for adaptive radiation. 

Additional files 



Additional file 1: References of morphological data used in this 
study. 

Additional file 2: GenBank accession number of the sequences 
used in this study. 

Additional file 3: Substitution model choice. 

Additional file 4: Phylogeny of the damselfish with branch lengths 
given in expected number of substitutions per site. 

Additional file 5: Factorial map of the MCA analysis with 
eigenvectors for the axes 1 and 2. 

Additional file 6: Factorial map of the MCA analysis with 
eigenvectors for the axes 3 and 4. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

GL and NS conceived the study. GL, CS and RW performed the analyses and 
drafted the manuscript. NS supervised GL, coordinated the project and 
helped to draft the manuscript. PB and NZ supervised RW, participated in 
the coordination and helped to draft the manuscript. All authors read and 
approved the final manuscript. 

Acknowledgments 

This work was funded by the grant (CRSl 13-125240) from the Swiss National 
Science Foundation to N.S., P.B.P. and N.E.Z. C.A.S. received a scholarship 
from the Summer Undergraduate Research Programme at Lausanne 
University. We would like to thank A. Dubuis, N. Arrigo, L. Pellissier, A. 
Kostikova, P-A. Christin, M. Alfaro, N. Alvarez, L. Keller and two anonymous 
reviewers for their suggestions and comments on previous versions of the 
manuscript. We thank A. Grutter and K. Cheney for the pictures and 
Rattlebrained Designs for the clownfish drawing. This work received support 
from the Vital-IT facilities from the Swiss Institute of Bioinformatics. 

Author details 

^Department of Ecology and Evolution, Biophore, University of Lausanne, 
1015, Lausanne, Switzerland. ^Swiss Institute of Bioinformatics, Genopode, 
Quartier Sorge, 1015, Lausanne, Switzerland. ^Ecological and Evolutionary 
Genetics Lab / Coral Reef Ecology Lab, School of Biological Sciences, The 
University of Queensland, 4072, St Lucia, Australia. ^Landscape Dynamics, 
Swiss Federal Research Institute WSL, 8903, Birmensdorf, Switzerland. 



Litsios et at. BMC Evolutionary Biology 201 2, 12:21 2 
http://www.bionnedcentral.conn/1471 -21 48/1 2/212 



Page 14 of 15 



Received: 6 June 2012 Accepted: 27 October 2012 
Publislied: 2 Novennber 2012 



References 

1. Grant PR: Ecology and evolution of Darwin's finclies. Princeton (NJ): Princeton 
University Press; 1986:1-458. 

2. Givnish TJ: Adaptive radiations and molecular systematics: issues and 
approaches. In Molecular Evolution and Adaptive Radiation. Edited by 
Givnish TJ, Systma KJ. Cambridge, U.K: Cambridge Univ. Press; 1997:1-54. 

3. Losos JB: Adaptive Radiation, Ecological Opportunity, and Evolutionary 
Determinism American Society of Naturalists E. 0. Wilson Award Address 
* Am Nat 2010, 175:623-639. 

4. Schluter D: The Ecology of Adaptive Radiation. Oxford, U.K.: Oxford University 
Press; 2000. 

5. Losos JB, Mahler DL: Adaptive Radiation: The Interaction of Ecological 
Opportunity, Adaptation, and Speciation. In Evolution since Darwin: The 
First 150 Years. Edited by Bell MA, Futuyma DJ, Eanes WF, Levinton JS. 
Sunderland, MA: Sinauer As; 2010:381-420. 

6. Salzburger W: The interaction of sexually and naturally selected traits in 
the adaptive radiations of cichlid fishes. Mol Ecol 2009, 18:169-185. 

7. Erwin DH: Lessons from the past: biotic recoveries from mass extinctions. 
Proc Natl Acad Sci USA 2001 , 98:5399-5403. 

8. Hughes C, Eastwood R: Island radiation on a continental scale: 
exceptional rates of plant diversification after uplift of the Andes. Proc 
Natl Acad Sci USA 2006, 103:10334-10339. 

9. Yoder JB, Clancey E, DES Roches S, Eastman JM, Gentry L, Godsoe W, Hagey 
TJ, Jochimsen D, Oswald BP, Robertson J, Sarver BJ, Schenk JJ, Spear SF, 
Harmon LJ: Ecological opportunity and the origin of adaptive radiations. 
JEvol Biol 2010, 23:1581-1596. 

10. Matschiner M, Hanel R, Salzburger W: On the origin and trigger of the 
notothenioid adaptive radiation. PLoS One 201 1, 6:e1891 1. 

11. Janson EM, Stireman JO, Singer MS, Abbot P: Phytophagous insect- 
microbe mutualisms and adaptive evolutionary diversification. Evolution 

2008, 62:997-1012. 

12. Gomez JM, Verdu M: Mutualism with Plants Drives Primate Diversification. 

Syst Biol 2012,61:567-577. 

13. Munday PL, van Herwerden L, Dudgeon CL: Evidence for sympatric 
speciation by host shift in the sea. Curr Biol 2004, 14:1498-1504. 

14. Bolnick Dl, Fitzpatrick BM: Sympatric speciation: models and empirical 
evidence. Annu Rev Ecol Evol Syst 2007, 38:459-487. 

15. Yoder JB, Nuismer SL: When does coevolution promote diversification? 
Am. Nat 2010, 176:802-817. 

16. Lengyel S, Gove AD, Latimer AM, Majer JD, Dunn RR: Ants sow the 
seeds of global diversification in flowering plants. PLoS One 2009, 
4:e5480. 

17. Quenouille B, Bermingham E, Planes S: Molecular systematics of the 
damselfishes (Teleostei: Pomacentridae): Bayesian phylogenetic analyses 
of mitochondrial and nuclear DNA sequences. Mol Phylogenet Evol 2004, 

31:66-88. 

18. Elliott JK, Mariscal RN: Coexistence of nine anemonefish species: 
differential host and habitat utilization, size and recruitment. Mar Biol 

2001, 138:23-36. 

19. Fautin DG, Allen GR: Anennonefishes and their host sea anennones. Perth: 
Western Australian Museum; 1997. 

20. Ollerton J, McCollin D, Fautin DG, Allen GR: Finding NEMO: nestedness 
engendered by mutualistic organization in anemonefish and their hosts. 
Proc R Soc B 2007, 274:591-598. 

21. Ricciardi F, Boyer M, Ollerton J: Assemblage and interaction structure of 
the anemonefish-anemone mutualism across the Manado region of 
Sulawesi. Indonesia Environ Biol Fishes 2010, 87:333-347. 

22. Mebs D: Chemical biology of the mutualistic relationships of sea 
anemones with fish and crustaceans. Toxicon 2009, 54:1071-1074. 

23. Buston PM, Garcia MB: An extraordinary life span estimate for the clown 
anemonefish Amphiprion percula. J Fish Biol 2007, 70:1710-1719. 

24. Fautin DG: The clownfish sea anemones: Stichodactylidae (Coelenterata: 
Actinaria) and others sea anemones symbiotic with Pomacentrid fishes. 
TAm Philos Soc 1981, 71:3-115. 

25. Puebia 0: Ecological speciation in marine v. freshwater fishes. J Fish Biol 

2009, 75:960-996. 

26. Jones GP, Planes S, Thorrold SR: Coral reef fish larvae settle close to home. 

Curr Biol 2005, 15:1314-1318. 



27. Almany GR, Berumen ML, Thorrold SR, Planes S, Jones GP: Local 
replenishment of coral reef fish populations in a marine reserve. Science 
2007,316:742-744. 

28. Parmentier E, Colleye 0, Fine ML, Frederich B, Vandewalle P, Herrel A: 
Sound production in the clownfish Amphiprion clarkii. Science 2007, 
316:1006. 

29. Colleye 0, Vandewalle P, Lanterbecq D, Lecchini D, Parmentier E: 
Interspecific variation of calls in clownfishes: degree of similarity in 
closely related species. BMC Evol Biol 201 1, 1 1:365. 

30. Parmentier E, Lagardere JP, Vandewalle P, Fine ML: Geographical Variation 
in Sound Production in the Anemonefish Amphiprion akallopisos. Proc R 
Soc B 2005, 272:1697-1703. 

31. Harmon LJ, Schulte JA, Larson A, Losos JB: Tempo and mode of 
evolutionary radiation in iguanian lizards. Science 2003, 301:961-964. 

32. Allen GR, Drew J, Kaufman L: Amphiprion barberi, a new species of 
anemonefish (Pomacentridae) from Fij'i, Tonga, and Samoa. Aqua 2008, 
14:105-114. 

33. Allen GR: Damselfish of the world. Melle, Germany: Mergus Publishers; 
1991:271. 

34. Elliott JK, Lougheed SC, Bateman B, Mcphee LK, Boag PT: Molecular 
evolution for the phylogenetic evidence of specialization in 
anemonefishes. Proc R Soc B 1999, 266:677-685. 

35. Frederich B, Fabri G, Lepoint G, Vandewalle P, Parmentier E: Trophic niches 
of thirteen damselfishes (Pomacentridae) at the Grand Recif of Toliara. 
Madagascar Ichthyol Res 2009, 56:1 0-1 7. 

36. Clabaut C, Bunje PME, Salzburger W, Meyer A: Geometric morphometric 
analyses provide evidence for the adaptive character of the Tanganyikan 
cichlid fish radiations. Evolution 2007, 61:560-578. 

37. Langerhans RB, Layman CA, Langerhans AK, Dewitt TJ: Habitat-associated 
morphological divergence in two Neotropical fish species. Biol J Linn Soc 
2003, 80:689-698. 

38. Drucker EG, Lauder GV: Locomotor function of the dorsal fin in teleost 
fishes: experimental analysis of wake forces in sunfish. J Exp Biol 2001, 
204:2943-2958. 

39. Kahilainen KK, Siwertsson A, Gjelland K0, Knudsen R, Bohn T, Amundsen P- 
A: The role of gill raker number variability in adaptive radiation of 
coregonid fish. Evol Ecol 2010, 25:573-588. 

40. Allen GR: Damselfishes of the south seas. Neptune city, N.J: T.F.H. 
Publications; 1975:240. 

41. Cowman PF, Bellwood DR: Coral reefs as drivers of cladogenesis: 
expanding coral reefs, cryptic extinction events, and the development of 
biodiversity hotspots. J Evol Biol 201 1, 24:2543-2562. 

42. Litsios G, Pellissier L, Forest F, Lexer C, Pearman PB, Zimmermann NE, 
Salamin N: Trophic specialization influences the rate of environmental 
niche evolution in damselfishes (Pomacentridae). Proc R Soc B 2012, 
279:3662-3669. 

43. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and 
high throughput. Nucleic Acids Res 2004, 32:1792-1797. 

44. Talavera G, Castresana J: Improvement of phylogenies after removing 
divergent and ambiguously aligned blocks from protein sequence 
alignments. Syst Biol 2007, 56:564-577. 

45. Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by 
sampling trees. BMC Evol Biol 2007, 7:214. 

46. Paradis E, Claude J, Strimmer K: APE: Analyses of Phylogenetics and 
Evolution in R language. Bioinformatics 2004, 20:289-290. 

47. R Development Core Team: R: A Language and Environment for 
Statistical Computing. Vienna, Austria: R Foundation for Statistical 
Computing; 201 1. 

48. Bellwood DR, Sorbini L: A review of the fossil record of the 
Pomacentridae (Teleostei: Labroidei) with a description of a new genus 
and species from the Eocene of Monte Bolca. Italy Zool J Linn Soc 1 996, 
117:159-174. 

49. Santini F, Harmon LJ, Carnevale G, Alfaro ME: Did genome duplication 
drive the origin of teleosts? A comparative study of diversification in ray- 
finned fishes. BMC Evol Biol 2009, 9:1 94. 

50. Near TJ, Eytan Rl, Dornburg A, Kuhn KL, Moore Ja, Davis MP, 
Wainwright PC, Friedman M, Smith WL: Resolution of ray-finned fish 
phylogeny and timing of diversification. Proc Natl Acad Sci USA 
2012, 109:13698-13703. 

51. Litsios G, Salamin N: Effects of Phylogenetic Signal on Ancestral State 
Reconstruction. Syst Biol 2012, 61:533-538. 



Litsios et at. BMC Evolutionary Biology 201 2, 12:21 2 
http://www.bionnedcentral.conn/1471 -21 48/1 2/212 



Page 15 of 15 



52. FitzJohn RG, Maddison WP, Otto SP: Estimating trait-dependent speciation 
and extinction rates from incompletely resolved phylogenies, Syst Biol 

2009, 58:595-611. 

53. Maddison WP, Midford PE, Otto SP: Estimating a binary character's effect 
on speciation and extinction. Syst Biol 2007, 56:701-710. 

54. Blomberg SP, Garland T, Ives AR: Testing for phylogenetic signal in 
comparative data: behavioral traits are more labile. Evolution 2003, 57:717-745. 

55. Pagel M: Inferring the historical patterns of biological evolution. Nature 
1999, 401:877-884. 

56. Revel I LJ: phytools: an R package for phylogenetic comparative biology 
(and other things). Methods Ecol Evol 201 2, 3:21 7-223. 

57. Orme CDL, Freckleton RP, Thomas GH, Petzoldt T, Fritz S: coper: Comparative 
Analyses of Phylogenetics and Evolution in R; 201 1. http://R-Forge.R-project. 
org/projects/caper/. 

58. Freckleton RP, Harvey PH: Detecting non-Brownian trait evolution in 
adaptive radiations. PLoS Biol 2006, 4:e373. 

59. Agrawal AA, Fishbein M, Halitschke R, Hastings AP, Rabosky DL, Rasmann S: 
Evidence for adaptive radiation from a phylogenetic study of plant 
defenses. Proc Natl Acad Sci USA 2009, 106:18067-18072. 

60. O'Meara BC, Ane C, Sanderson MJ, Wainwright PC: Testing for different 
rates of continuous trait evolution using likelihood. Evolution 2006, 
60:922-933. 

61 . Cooper JW, Smith LL, Westneat MW: Exploring the radiation of a diverse reef 
fish family: phylogenetics of the damselfishes (Pomacentridae), with new 
classifications based on molecular analyses of all genera. Mol Phylogenet Evol 

2009, 52:1-16. 

62. Marko PB: Fossil calibration of molecular clocks and the divergence times 
of geminate species pairs separated by the isthmus of panama. Mol Biol 
Evol 2002, 19:2005-2021. 

63. Lessios HA: The Great American Schism: Divergence of Marine Organisms 
After the Rise of the Central American Isthmus. Annu Rev Ecol Evol Syst 

2008, 39:63-91. 

64. Alfaro ME, Santini F, Brock C, Alamillo H, Dornburg A, Rabosky DL, Carnevale 
G, Harmon LJ: Nine exceptional radiations plus high turnover explain 
species diversity in jawed vertebrates. Proc Natl Acad Sci USA 2009, 
106:13410-13414. 

65. Nee S: Mooers a 0, Harvey PH: Tempo and mode of evolution revealed 
from molecular phylogenies. Proc Natl Acad Sci USA 1 992, 89:8322-8326. 

66. Rabosky DL: Extinction rates should not be estimated from molecular 
phylogenies. Evolution 2010, 64:1816-2184. 

67. Revell LJ, Harmon LJ, Collar DC: Phylogenetic signal, evolutionary process, 
and rate. Syst Biol 2008, 57:591-601. 

68. Ackerly D: Conservatism and diversification of plant functional traits: 
Evolutionary rates versus phylogenetic signal. Proc Natl Acad Sci USA 

2009, 106(Suppl):1 9699-1 9706. 

69. Losos JB: Lizards in an evolutionary tree: ecology and adaptive radiation of 
anoles. Berkeley: University of California Press; 2009:528. 

70. Bellwood DR, Wainwright PC: Locomotion in labrid fishes: implications for 
habitat use and cross-shelf biogeography on the Great Barrier Reef. Coral 
/?eefs 2001, 20:139-150. 

71. Lovette IJ, Bermingham E, Ricklefs RE: Clade-specific morphological 
diversification and adaptive radiation in Hawaiian songbirds. Proc R Soc B 
2002, 269:37-42. 

72. Harmon LJ, Melville J, Larson A, Losos JB: The role of geography and 
ecological opportunity in the diversification of day geckos (Phelsuma). 
Syst Biol 2008, 57:562-573. 

73. Pushparaj A: Intra specific hybridization between Amhiprion sebae 
and A. polymnus under captive conditions. Int J Biological Technology 

2010, 1:52-56. 

74. Thresher RE, Colin PL, Bell U: Planktonic Duration, Distribution and Population 
Structure of Western and Central Pacific Damselfishes (Pomacentridae). 

Copeia 1989, 1989:420-434. 



doi:1 0.1 186/1471-2148-12-212 

Cite this article as: Litsios et al.: Mutualism with sea anemones triggered 
the adaptive radiation of clownfishes. BMC Evolutionary Biology 2012 
12:212. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



